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Abstract 

Properties of the g-ultraspherical polynomials for q being a primitive root of unity are derived using 
a formalism of the S0g(3) algebra. The orthogonality condition for these polynomials provides a new 
class of trigonometric identities representing discrete finite-dimensional analogs of q-beta integrals of 
^ ' Ramanujan. 

m 

m . Mathematics Subject Classifications (1991). 17B37, 33D80 

I 1. Introduction. There are many applications of quantum algebras, or g-analogues of Lie algebras in 

' mathematical physics Q. In particular, it is well known that the corresponding representation theory is 
■ related to g-special functions and g-orthogonal polynomials. In this Letter we investigate from this point 
of view a class of g-ultraspherical polynomials for q a root of unity. 

The monic g-ultraspherical polynomials satisfy three-term recurrence relation 

Pn+lix) + UnPn-lix) = xPn{x), 71 = 1,2,... (1) 

I and initial conditions Po{x) = 1, -Pi (a;) = x, where coefficients Un have the form IH, |^, ^ 

^ . (1 -g")(l 

(1 - g"+/3)(l - q 



. „„+/3wi „„+/3_n- l^J 



These polynomials are well investigated for < g < 1. The case of g a root of unity was considered 
only in the nontrivial limits /? — > or 1 leading to the so-called sieved ultraspherical polynomials j^. In 
these limits the recurrence coefficients Un > for all n > 0, so that the sieved polynomials are orthogonal 
on an interval of real line with some continuous measure. Consider what happens if /? 7^ 0, 1 and q is 
A^th primitive root of unity, q = exp 2'iTip/N, (p, N) = 1. Then ti„ can be rewritten in the form 

sincjn sina;(n + 2/3 — 1) 
— (3) 



sinaj(n + f3) sina;(n + /? — 1) ' 

where u = irp/N. In order for to possess a positive weight function (i.e. to be "classical" orthogonal 

polynomials) ti„ must be positive. Reality of Un implies reality of the parameter p. Then (^) cannot 

^Published in Lett. Math. Phys. 37 (1996), pp. 173-180. 
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satisfy the positivity condition for all n > 0. One has to analyze thus the finite-dimensional cases when 
uq = um = for some natural M > 2 and > for < n < M (then M is the dimension of Jacobi 
matrix defining the eigenvalue problem (|l|)). 

The possibility to construct finite-dimensional classical orthogonal polynomials on trigonometric grids 
(including those generated by g a root of unity) was mentioned in however no explicit examples 
were presented. For examples of non-standard (/-special functions at = 1 related to the Lame-type 
differential equations, see 0. 

Recently we have found that a special class of g-ultraspherical polynomials for = 1 can be obtained 
by "undressing" of the finite-dimensional Chebyshev polynomials |9| , existence of the more complicated 
cases has been pointed out as well. Here we describe g-ultraspherical polynomials with positive measure 
when q is the p = 1 primitive root of unity, i.e. when 

q = exp27ri/N, iV = 2,3, .... (4) 

The main tool in our analysis will be the theory of representations of quantum algebra S0q{3) intro- 



duced in [1^, 11 1 . The connection of this algebra with some (finite-dimensional) class of g-ultraspherical 
polynomials for real q was mentioned in |l^ . 

2. Restrictions upon the Parameters. Consider restrictions upon the parameters M and /3 necessary 
for positivity of ti„ in the interval < n < M. There are two principally different situations. The first 
possibility is M = A^, when dimension of representation coincides with A^. Then the positivity condition 
leads to the following restrictions upon /3 (all the inequalities are modulo A^): 

-l/2</?<0, 0</3<l, l</?<3/2. (5) 

Note that the ends of intervals do not belong to this class. Indeed, there are ambiguities in choosing the 
coefficients Un for P = 0,1. The cases /? = —1/2, 3/2 lead to um = for M = 2, N — 2, which is not 
allowed since we demand the dimension of representation to be equal to A^. 
The second possibility arises if 2/3 is a natural number: 

2/3 = j = l,2,3,...,A^-l. (6) 

We define coefficients Un in the case /3 = 1 as ?/o = un-i =0, ui = . . . = um-2 = 1- Then the dimension 
of representations defined by (^) is determined by the general formula M = A^-|- 1 — j. We note by passing 



that the even j cases describe undressing of the "discrete finite well" in discrete quantum mechanics |13]. 
The case j = 1 corresponds to the (/-Legendre polynomials because Un reduce to the ordinary Legendre 
polynomials' recurrence coefficients for N ^ oo. This system has an interesting physical application in 
the Azbel-Hofstadter problem of electron on a two-dimensional lattice in a magnetic field p^ . 

It is not difficult to see that (^) and (^ are the only choices of P providing positive Un (we skip the 
trivial case /? = —1/2 when M = 2 and assume /3 7^ 0). 

3. sOg(3) Algebra and Its Representations. Consider the associative algebra of three generators defined 
by commutation relations 

[Ko,Ki]^ = K2, [Ki,K2]^ = -Ko, [K2, Kq]^ = -K^, (7) 

where [A, B]^ = e^^/'^AB — e~''-^^'^BA denotes so-called (7- mutator and oj = tt/N. Note that if Kq and Ki 
are hermitian operators then K2 should be anti-hermitian. Note also that for N ^ 00 we get the ordinary 
rotation algebra so(3), which justifies the name S0g(3) for the algebra (0). This algebra was introduced (in 
a slightly different form) in 1986 by Odesskii |l^ who considered its representations without discussion of 



unitarity. Its relations to quantum algebra suq{2) and g-special functions have been discussed in |11, 12] 
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A non-compact version of this algebra S0q{2, 1) for real q appeared to be quite useful as an algebra 
describing dynamical symmetry of special classes of discrete reflectionless potentials 

Construct the unitary finite-dimensional representations of S0g(3) by the method, proposed in |12]. 
It is easy to see that there exists an orthonormal basis \n) for which the operator Kq is diagonal whereas 
the operator Ki is two-diagonal: 

Ko\n) = Xn\n) , = a„+i|n + 1) + a„|n - 1), (8) 

where A„ are eigenvalues of operator Kq and a„ are matrix elements of a representation. Substituting 
dD into (I) we find 

^^^ cosc.(n + /?) ^ n = 0,l,...,M-l, (9) 
smo; 

2 — cosa;(n + /3) cosa;(n + /3 — 1)) ^^^^ 

" 4sin^ w sina;(n + /3) sina'(n + /3 — 1) ' 

where v is the eigenvalue of Casimir operator of soq{3) algebra which can be represented in the form 

Q = ^{K2k2 + K2K2) - cosio{Ki + Kf), (11) 

where K2 = [i^Oj -^ij-i^;- The condition ao = gives = cos lj (3 cos — 1), and then 

2 sinwn sina;(n + 2/3 — 1) ^^^^ 

" 4sin^ti;sina;(n + /3) sinti;(n + /? — 1) 

coincide up to a constant factor with (|3|). Hence, all unitary irreducible finite-dimensional representations 
of the S0q{3) algebra for q = exp27ri/N are exhausted by two possibilities: @ of the dimension A'^ and 
(|6|) of a smaller dimension. 

We introduce for brevity the following terminology: the representations (^) will be called "comple- 
mentary series" ; the cases (|6|) will be called "integer series" for /? = 1, 2, 3, . . . and "half-integer series" for 
P = 1/2,3/2, . . . (the case /3 = 1/2 is exceptional, for it belongs to both complementary and half-integer 
series). Such definitions are supported by the fact that eigenvalues of the Casimir operator (0) are 
continuous and quantized for the complementary and integer (half-integer) series respectively (cf. with 
the Lorentz algebra). 

4. q- Ultraspherical Polynomials as Overlap Coefficients. g-Ultraspherical polynomials arise as overlap 
coefficients between two distinct bases in the space of S0q(3) representations. Indeed, because Ki is 
hermitian, we can choose another orthonormal basis |s) for which the operator Ki is diagonal: 

Ki\s) = fis\s). (13) 

sOq{3) is a specisal case of the AT^(3) algebra describing symmetries of the Askey- Wilson polynomials 
||T5[| . From general properties of the latter algebra we know that the operator Kq cannot be more than 
tridiagonal in this dual basis: 

Ko\s) = ds+i\s + 1) + 4|s - 1) + bs\s), (14) 

where s = 0, . . . ,M — 1. The explicit form of matrix elements dg and bg depends on the representation 
series of sOq(3). For integer and half-integer series we have the spectrum 

= ^°^"^^+^'/'\ (15) 
sinw 
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where j = 2/3 = 1, 2, 3, . . . , and the matrix elements 



4= , sin^.sin^.+j-l) ^^^^ ^^g^ 

4 sin a; sinu;(s + j'/2) sina;(,s — 1 + j'/2) 

This case is symmetric - matrix elements in the bases |n) and \s) are identical. 
For the complementary series we have 

^^^cosuAs + m^ . = 0,l,...,iV-l, (17) 



smtij 

sina;(s - /3 + 1/2) sinu;(s + /3 - 1/2) 
4sin^a;sina;(s + l/2)sina;(s - 1/2) ' 



s = l,2,...,iV-2, (18) 



^0 = ^N-l = ^^'^^^(^/^ = = hi=h2 = ... = bN-2 = 0. (19) 

2 smc<j sma;/2 

Note that the coefficients dn do not truncate automatically at the ends of the index intervals, one has to 
do it by hands. In this case the operator Kq is tridiagonal instead of being two-diagonal as it is expected 
from the symmetry of S0q{3) (after the replacement K2 ii^2 this algebra looks totally symmetric 
under the cyclic permutations). This anomaly shows that permutations of generators are not necessarily 
unitary automorphisms of the "Cartesian" quantum algebras. 

For any series: integer, half-integer or the complementary one, we can find the overlap coefficients 
between two bases {s\n). These coefficients can be factorized, (s|n) = (s|0)5„(/is), where SnitJ-s) are 
symmetric polynomials satisfying the three-term recurrence relation 

O-n+lSn+l + a,nSn-l = fJ-sSri, Sq = 1, Si = fJ-s/ai- (20) 

Polynomials 5„(/is), n > 1 are connected with Pn{x) (|l|) by the simple relation 

Pnix) = ^/UiU2 ■ ■ ■ UnSnif^s), X = 2flsSmuJ. (21) 

Existence of the algebraic interpretation of g-ultraspherical polynomials in terms of the sOq{3) algebra 
allows to calculate the weight function for these polynomials. Omitting the details (we use the method 
described, e.g. in |l^), we present the final result. 

For the integer and half-integer series the weight function has the form 

i-i 

WsU) =smio{s + j/2)Y[smuj{s + l), s = 0,1, . . . , N - j. (22) 

1=1 

For the complementary series we have 

, sina;(s + l/2) Vt^ sina;(/? + 1/2 + /) , , 

Ws{2P)=wo ' ' U- / /Q^o/o^ V s = l,2,...,N-l, 23 

sma;/2 ^-^-^ sma;(— p -|- 3/2 + <) 

where wq is (undefined) value of the weight function at the points ,5 = and s = N — 1. 

Clearly both weight functions (^) and ( p3|) are known up to a normalization constant. This constant 
can not be found directly from the representation theory of sOq{3) algebra and should be calculated 
separately. For the integer and half-integer series it is possible to derive weight functions together with 
normalization constants with the help of Darboux transformations. 
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5. Darboux Transformation and Normalization Constants. Let Pn\x) be the monic g-ultraspherical 
polynomials belonging to integer or half-integer series. Write the orthogonality relation in the form 

N-j 

J2 P^'\xs)P^Hxs)Ws{j) = K{j)5nm, (24) 
s=0 

where hn{j) are normalization constants which have to be found, and Ws{j) is given by (^). Recall that 
for the taken series Xs{j) = 2cosu{s + j/2). 

The crucial observation is that the polynomials Pn^'^\x) and Pn\x) are related to each other by 
the Darboux transformation which is equivalent in our case to the transition to symmetric kernel 
polynomials [16|: 

+ 2)) = (25) 

where 

AM ^ (26) 



From the theory of kernel polynomials |16| | it follows that the weight function is transformed as 

ws{j + 2) = ws+i{j){xl{j) - x2^i(j))/4, (27) 
and the normalization constants are transformed as 

hnij + 2) = Kij)An{j)/A. (28) 

In our case the factor A„(j) can be easily found 

Anij) = (29) 

sina;(n + 1 + j/2) sinti;(ra + j/2) 



Then using (29), ( p^ ) and starting from j = 2 and j = 1 we obtain explicit expressions for the normal- 
ization constant in the case of integer series, j = 2k: 

h^{2k) = ^'0{2)jn+k+lSn+k+2 ■ • • gn+2fe-l ^ ^ = 0, 1 , . . . , iV - 2A;, (30) 
Sn+lSn+2 ■ ■ ■ Sn+k-1 

where s„ = sinwn, and in the case of half-integer series, j = 2k + 1: 

Ki2k + 1) = 2 ' - • • , n = 0,l,...,iV-2A:-l, (31) 

4 Sl/2S3/2'S5/2 • • • ■5„+fc_l/2'5n+fc+l/2 

where we assume that for n = the product ■ ■ ■ is replaced by 1. The only coefficients left to be 
determined are /io(2) and /io(l)- But these constants can be calculated directly. Indeed, for j = 2 we 
have the finite-dimensional Chebyshev polynomials 

PniXs) = -^"--(-+^)(;+^) . (32) 

sm L(j(s -|- Ij 

The weight function in this case is Ws{2) = sin^ u>{s + 1). Hence 

N-2 

ho{2) = J2 sin^ ^(« + 1) = ^/2- (33) 

s=0 
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The g-Legendre polynomials case j = 1 is characterized by the surprisingly simple weight function, 
Ws{l) = sinu;(s + 1/2), from which we derive 



N-l 

ho{l) = E smu;(5 + 1/2) = ——. (34) 

s=0 ' 



Thus we have calculated hn{j) for both integer (|3^) and half-integer (|3l|) series. 

The n = values of the normalization constants hn{j) provide non-trivial trigonometric identities of 
the form X^^o^ Wr{j) = h^ij)- For j = 2k we get 



N-2k riAT 

^Jv Sfc+iSfc+2 • • • S2k-1 



Sr+kSr+lSr+2 ■ ■ ■ Sr+2k-l = TT , i-jt)J 

4'=SiS2...Sfe-l 

where A; = 2, 3, ... , [A^/2]. For j = 2A; + 1 we get 

N-2k-l 

E_ SiS2 ■ ■ ■S2k /„„x 
Sr+k+l/2Sr+lSr+2 ■ ■ ■ Sr+2k — 2 2 ' v"^"/ 

r=0 ^ '^1/2^3/2 • • • '^fc-l/2'^fc+l/2 

where = 1,2, ... , [A^/2]. These identities can be considered as discrete finite-dimensional analogs of 
the g-beta integrals of Ramanujan (the latter are defined for real q, see e.g. p). Of course, in the 
case of (7-ultraspherical polynomials we get only special cases of such identities. Their generalizations are 
obtained if we deal with general trigonometric (four-parameter) analogs of the Askey- Wilson polynomials. 
Existence of similar simple relations for the representations of complementary series is an interesting open 
problem which we hope to address in the future. 
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